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Abstract 

A new model is proposed to a collapsing star consisting of an anisotropic fluid with bulk viscosity, 



5-H 

radial heat flow and outgoing radiation. In a previous paper one of us has introduced a function 
time dependent into the g rr , besides the time dependent metric functions ggo and g^. The aim 
of this work is to generalize this previous model by introducing bulk viscosity and compare it to 

adiabatic index is analyzed. Our work is also compared to the case of a collapsing fluid with 

o 

bulk viscosity of another previous model, for a star with 6 M . The pressure of the star, at the 
beginning of the collapse, is isotropic but due to the presence of the bulk viscosity the pressure 



^ , becomes more and more anisotropic. The black hole is never formed because the apparent horizon 

H ■ 

formation condition is never satisfied, in contrast of the previous model where a black hole is 
formed. An observer at infinity sees a radial point source radiating exponentially until reaches 
the time of maximum luminosity and suddenly the star turns off. In contrast of the former model 
where the luminosity also increases exponentially, reaching a maximum and after it decreases until 
the formation of the black hole. The effective adiabatic index diminishes due to the bulk viscosity, 
thus increasing the instability of the system, in both models, in the former paper and in this work. 
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I. INTRODUCTION 



One of the most outstanding problems in gravitation theory is the evolution of a col- 
lapsing massive star, after it has exhausted its nuclear fuel. The problem of constructing 
physically realistic models for radiating collapsing stars is one of the aims of the relativistic 
astrophysics. However, in order to obtain realistic models we need to solve complicated 
systems of nonlinear differential equations. In many cases we can simplify the problem 
considering some restrictions in these equations and solve the system analytically. Such 
models, although simplified, are useful to construct simple exact models, which are at least 
not physically unreasonable. This allows a clearer analysis of the main physical effects at 
play, and it can be very useful for checking of numerical procedures. 

A great number of the previous works in gravitational collapse have considered only 
shear-free motion of the fluid [l-4|. This simplification allows us to obtain exact solutions 
of the Einstein's equations in some cases but it is somewhat unrealistic. It is also unrealistic 
to consider heat flow without viscosity but if viscosity is introduced, it is desirable to allow 
shear in the fluid motion. 

In the work [5;] the authors have studied the collapse of a radiating star with bulk vis- 
cosity but they still maintained the shear-free motion of the fluid. Thus, it is interesting to 



study solutions that contains shear,because it 
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P 



12 



ays a very important role in the study of 



13) and in 



gravitational collapse, as shown in p, ]7|, ]9|, 

In the first paper jfj, [7] we have compared two collapsing model: a shear- free and a 
shearing model. In this model we have imposed that the metric components g u and g rr were 
independent of the time and only g ee and were time dependent. We were interested in 
studying the effect of the shearing motion in the evolution of the collapse. It was shown 
that the pressure of the star, at the beginning of the collapse, is isotropic but due to the 
presence of the shear the pressure becomes more and more anisotropic. The anisotropy in 
self-gravitating systems has been reviewed and discussed the causes for its appearance by 



Herrera and Santos 151 ] . As shown by Chan 



the simplest cause of the presence of 



anisotropy in a self-gravitating body is the shearing motion of the fluid, because it appears 
without an imposition ad- hoc [41. 

q n q 

In the second work |9[ we have used the same model of Chan [6|, |7| and we have analyzed 
a collapsing radiating star consisting of an anisotropic fluid with shear viscosity undergoing 
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radial heat flow with outgoing radiation, but without bulk viscosity. 

In the third paper [h| we have also used the same model previous papers j^] j3] and 
we have analyzed a collapsing radiating star consisting of an anisotropic fluid with bulk 
viscosity undergoing radial heat flow with outgoing radiation, but without shear viscosity 

In the fourth work ll| we have generalized our previous models by introducing a function 



time dependent into the g rr and in a recent paper 13J we have introduced the shear viscosity. 



The aim of this work is to generalize our previous model [10] by introducing a time 
dependent function into the g rr , besides the time dependent metric functions gee and g^, 
and to compare the physical results with the previous ones. 

This work is organized as follows. In Section 2 we present the Einstein's field equations. 
In Section 3 we rederive the junction conditions, since in the former paper [l3| have obtained 
only results without bulk viscosity. In Section 4 we present the proposed solution of the 
field equations. In Section 5 we describe the model considered in this work for the initial 
configuration. In Section 6 we present the energy conditions for a bulk viscous anisotropic 
fluid. In Section 7 we show the time evolution of the total mass, luminosity and the effective 
adiabatic index and in Section 8 we summarize the main results obtained in this work. 



II. FIELD EQUATIONS 

We assume a spherically symmetric distribution of fluid undergoing dissipation in the 
form of heat flow. While the dissipative fluid collapses it produces radiation. The interior 
spacetime is described by the most general spherically symmetric metric, using comoving 
coordinates, 



ds 2 



-A 2 (r, t)dt 2 + B 2 {r, t)dr 2 + C 2 (r, t){d9 2 + sin 2 



(1) 



The exterior spacetime is described by Vaidya's 16[ metric, which represents an outgoing 
radial flux of radiation, 



ds z 



2m(v) 



dv - 2dvdr + r 2 (d6 + sin 



(2) 



where m(v) represents the mass of the system inside the boundary surface E, function of 
the retarded time v. 
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We assume the interior energy-momentum tensor is given by 



G a/3 = KT a/3 = k[(h + pt)u a up + p t g a i3 + (p~ Pt)X a Xp 

+q a up + q p u a - C© (g a p + u a up)} , (3) 

where fj, is the energy density of the fluid, p is the radial pressure, p t is the tangential 
pressure, q a is the radial heat flux, X a is an unit four-vector along the radial direction, u a 
is the four-velocity, which have to satisfy u a q a = 0, X a X a = 1, X a u a = and k = 8n (i.e., 
c = G = 1). The quantity ( > is the coefficient of bulk viscosity and the shearing tensor 
a a p is defined as 

<7a/3 = «(a;/9) + «(a«/3) - \^>{9ap + U a Up), (4) 

with 

u a = u a . p u p , (5) 

© = < Q , (6) 

where the semicolon denotes a covariant derivative and the parentheses in the indices mean 
symmetrizations. 

Since we utilize comoving coordinates we have, 

u a = A'X, (7) 

and since the heat flux is radial 

q a = q5«. (8) 
Thus the non-zero components of the shearing tensor are given by 

IB 2 (B C\ 

C 2 ( B C\ 

= °22 sin 2 6. (11) 

A simple calculation shows that 



Thus, we define the scalar a as 



a = — 



1 (B C s 
3A \B ~ C 



Using dTJ and (jSJ), we can write that 



The non- vanishing components of the field equations, using ([[]). 
interior of the boundary surface £ are 



(13) 



(14) 

©, (ED, © and dUD , 



G, 



oo 



A\ 2 



C \C C B 



+ 



A\ 2 (7/(7 B' 
CJ + C\C + B J 



M 2 /i, 



(15) 



a) 



2 C - + (^ 



2^ 

A c 



KJ B 2 ( P -Ce), 



(16) 



G 



22 





"C" 








2 











A" CM' A' 5' 

+ T + cl ~ Ib 

B~C~ CB 



KC 2 ( Pt -ce), 



B' C 

~b7j 



AC AB 
AC + AB 



(17) 



G 33 = G 22 sm 2 9, 



(18) 



G 



01 



(19) 



The dot and the prime stand for differentiation with respect to t and r, respectively. 
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III. JUNCTION CONDITIONS 



We consider a spherical surface with its motion described by a time-like three-space E, 
which divides spacetimes into interior and exterior manifolds. For the junction conditions 



we follow the approach given by Israel 
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181 ] . Hence we have to demand 



(ds 2 



(20) 



where is the extrinsic curvature to S, given by 



(21) 



, = , d 2 x a _ , drfdtf_ 



(22) 



and where are the Christoffel symbols, the unit normal vectors to E, x a are the 
coordinates of interior and exterior spacetimes and are the coordinates that define the 
surface E. 

From the junction condition ( 1201) we obtain 



(23) 



C{rv,t)=Tx{v), 



(24) 



' dv 
dr 



-2 



2m ^dr 
r dv 



where r is a time coordinate defined only on E. 



The unit normal vectors to E (for details see 19] ) are given by 



(25) 



(26) 



n. 



1 - 



2m 



2- 

dv , 



-1/2 



(27) 



The non-vanishing extrinsic curvature are given by 



6 



k: 



' dtyA'A 

dr 



B 



K a 



'C'C* 



K A , = K aa sin 2 8, 



(28) 
(29) 
(30) 



k: 



d 2 v ( dv 
dr 2 \ dr 



. dr i 



-i 



2m 
r 



dv\m 
dr / r 2 



dr 
r + — r 

dr 



(31) 



(32) 



K+=K+ sin 2 9. 



From the equations (1291) and (1321) we have 



(33) 



dr ) \ r / dr 



s \ / s 

With the help of equations ()23D, flSD, (ESJ), we can write flM]) as 



C 

"' '2 



which is the total energy entrapped inside the surface £ 
From the equations (128]) and (13T1) . using (|23|) . we have 



20). 



d 2 v ( dv^ 
d^Xd^ 



dv\m 
dr I v 2 



' AT 
AB 



Substituting equations ( l23l) . ( 1241) and ( 13~5l) into ( 1341) we can write 



dv 
dr 



'C 

B + I 



-i 



(34) 



(35) 



(36) 



(37) 

Differentiating ( 1371) with respect to r and using equations ( 1351) . ( 1371) . we can rewrite (1361) 



as 
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Comparing (1381) with 0161) and 0191) . we can finally write 



(39) 



This result is analogous to the one obtained by Chan 
now we have generalized for an interior fluid with bulk viscosity. 
The total luminosity for an observer at rest at infinity is 



ll| for a shearing fluid motion but 



dm\ 
dv /, 



dm dt f dv 
dt drXdr, 



(40) 



Differentiating fl35|) with respect to t, using fl23|) . f )37|) and f|T6|) . we obtain that 



/v 



(p - ce)c 2 



'a c" 
+ 



(41) 



B A) 

The boundary redshift can be used to determine the time of formation of the horizon. 
The boundary redshift zs is given by 



dv 
dr 



1 + 



(42) 



The redshift, for an observer at rest at infinity diverges at the time of formation of the 
black hole. From fl37|) we can see that this happens when 



B 



A 



0. 



(43) 
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IV. SOLUTION OF THE FIELD EQUATIONS 



Again as in former paper we propose solutions of the field equations (|T5|) -( IT9|) with 
the form 

A(r,t) = A (r), (44) 



B{r,t) = B (r)h(t), 



(45) 



C(r,t) = rB (r)f(t). 



(46) 



We have chosen this separation of variables in the metric functions, in order to have the 
following properties: (a) when h(t) — > 1 and f(t) — > 1 the metric functions represent the 
static solution of the initial star configuration; (b) when h(t) = f(t) the metric functions 
represent the shear-free solution. We also remark that, following the junction condition 
equation (|24|) . the function C(rs,t) represents the luminosity radius of the body as seen by 
an exterior observer. On the other hand, with this solution the proper radius f r B(r\t)dr 
evolve with the time. Such a property was not present in the previous models {?! |Q[| flO| 
12|- 



Thus, the expansion scalar f)T4j) can be written as 

9 

Now the equations (|T5|) - (|T9|) can be written as 



(47) 



^ + M0(/ +2 ^) + ^o 2 (/ 12 h\ 



'h 2 



(48) 



Po 1 



1 (\ 1 \ < (2f 
r 2 Bl\P h 2 ) A \f + h 



(49) 



Po 



i (f 



hf 



h 2 



Al \f + h + hf / 



Ao\f h J 



(50) 
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/■•'/ 



fh\ , 



r 



r 



1 

r 



(51) 



where 



Wo = -752 



1 

m 



'Bo 



B rB 



npo 



R2 
-°0 



B n 



2 -£?n -^4.n Bc\ 2 ./In 

r Bo Ao Bo r A 



(52) 



(53) 



We can see from equations (|4"8"|) - (|5"T|) that when the functions h(t) = 1 and f(t) = 1 we 
obtain the static perfect fluid configuration. 

Substituting equations (|4"9~j) and ( 15"TT) into ( 1391) . assuming also that po( r s) = 0, we obtain 
a second order differential equation in h(t) and /(£), 



2 HK 



B (/)-"(£ 



+ 6 



,/ 2 



/i 2 



(54) 



where 



5 j U r A , 



21^ f^ + i 



Br 



and 



5 n 



(55) 
(56) 

(57) 



In order to obtain the quantities (l48 i) -( l5Ti) we have first to find an appropriate h(t) 
function. Let us first assume that f(t) = 1. Thus, we obtain the differential equation for 
h(t) given by 



a h - h 2 + 1 = 0, 
where a = a/b and whose solution is given by [8] 



tanh I 1 

,a 



(58) 



(59) 
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Now, using equation (|58|) we can write equation f|54|) in the following way 



2ff + f 2 + a(f/h)f + b(l-f 2 ) = 0. (60) 

This equation is almost identical to the one obtained previously Q Q jsj 10] 12], except 
the factor (1/h) in the third term. Thus, as before it has to be solved numerically, assuming 
that at t — > — oo represents the static configuration with f(t — > — oo) —> and fit — > 
— oo) — > 1. We also assume that fit — > 0) —> 0. This means that the luminosity radius 
C(rE, t) has the value r^Boir^,) at the beginning of the collapse and vanishing at the end of 
the evolution. Analogously, the proper radius has the value hit) f£ B Q ir)dr at the beginning 
of the collapse and also vanishing at the end of the collapse. 

V. MODEL OF THE INITIAL CONFIGURATION 

We consider that the system at the beginning of the collapse has a static configuration 
of a perfect fluid satisfying the Schwarzschild interior solution j2l| 

A ° ~ 2(l + ri)(l + r>)' (61) 



2R 

1 + r 

where 



B o = 7-^> (62) 



g{r) = 3(1 - r|)(l + r 2 ) - (1 + r|)(l - r 2 ), (63) 

and 

R = r nv ^±f. (64) 

and where rs is the radial coordinate relative to the physical initial radius of the star in 
comoving coordinates and mo is the initial mass of the system. Thus the static uniform 
energy density and static pressure are given by 



= (65) 
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FIG. 1: Time behavior of the functions f(t) and h(t) for the model with or without bulk viscosity. 
The time is in units of second and f(t) and h(t) are dimensionless. The symbols £ > mean that 
the plotted quantity is independent of C- 
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FIG. 2: The function hf/f as a function of the time. The time is in units of second. The symbols 
( > mean that the plotted quantity is independent of (. 
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FIG. 3: The expansion scalar profiles as a function of the time. The radial coordinates r and 
are in units of second. 
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FIG. 4: Density profiles for the model with or without bulk viscosity. The radial coordinates r and 
rs are in units of seconds and the density is in units of sec~ 2 . The symbols ( > mean that the 
plotted quantity is independent of £. 
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FIG. 5: Radial pressure profiles for four different values of ( . The radial coordinates r and rs are 
in units of seconds and the radial pressure is in units of sec -2 . 
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FIG. 6: Tangential pressure profiles for four different values of £ . The radial coordinates r and 
are in units of seconds and the tangential pressure pt is in units of sec -2 . 
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FIG. 7: The profiles for four different values of £ of the ratio between the radial and tangential 
pressures. The radial coordinates r and rg are in units of seconds; and the radial and tangential 
pressure^ and pt, are in units of sec -2 . 
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FIG. 8: Heat flux scalar profiles for the model with or without bulk viscosity. The radial coordinate 
r and rg are in units of seconds and the heat flux q is in units of sec -2 . The symbols ( > mean 
that the plotted quantity is independent of 
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6 (ry — r 2 ) 
R 2 g(r) 

We consider the initial configuration as due to a iron core of a presupernova with mo = 
6M Q , initial radial coordinate = 1.6 x 10 5 km, which correspond to 2.963 x 10~ 5 and 
5.337 x 10 _1 , respectively, in units of second. Thus, the physical radius rx-BofVs) = 25.742 
km, which gives a density of 1.675 x 10 14 g cm -3 22]. With these values we can solve 
numerically the differential equation (|60|) . We can see from (JSTjl . using ( 161~]) -( 164")) and this 
initial configuration, that [{B' /B + 1/r - A' /A )/A ]x < 0, (B' /B + l/r) E > 0, g < 
and by the fact that > then we conclude that / < 0. In figured] we can see the time 
evolution of the functions f(t) and h(t). 

In order to determine the time of formation of the horizon /bh, we use the equations (1431) . 
dSD-dll, dMD-dSnD and write 

^h hh = _ 2 4( 1 -4) 2 w _ 3 6Q6 x 1q3 (67) 
/ bh m (l + r|;) 4 

Using the numerical solution of /(t), and equation ( |671) . we can see from figure 
[2] that the horizon is never formed, because the function hf/f does not reach the value 
—3.606 x 10 3 . At the first sight this fact could be interpreted as the formation of a naked 
singularity. However, this is not the case as we will see below in the calculation of the total 
energy entrapped inside the hypersurface E. 

We will assume that ( is constant, but in general the bulk viscosity coefficient depends 
on the temperature and density of the fluid [23| . The dependence of the expansion scalar on 
the time and radial coordinate is shown in the figure |3j Hereinafter, the values of ( will be 
1.347 x 10 30 , 6.736 x 10 30 and 1.347 x 10 31 g cm" 1 s _1 , which correspond to values 100, 500 
and 1000 s" 1 , respectively, in time units. These values are about ten orders of magnitude 
above current estimates of the bulk viscosity coefficient in neutron stars J24]. If these lower 
values were used in our model, we would have obtained results like the ones with ( m 0, i.e., 
without bulk viscosity. 

It is shown in figures H] and El the radial profiles of the density and the heat flux. It is 
shown only one plot for each quantity because they do not depend on the bulk viscosity, 
which can be seen from equations ( 148|) and (I5ip . 

In figure |5] and [6] we notice that the radial and tangential pressures diminish with the 
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bulk viscosity. 

In the figure [7] (C = 0) we can see that the star is isotropic at the beginning of the collapse 
(/ = 1) but becoming more and more anisotropic at later times. 



VI. ENERGY CONDITIONS FOR A VISCOUS ANISOTROPIC FLUID 



Following the same procedure used in Kolassis, Santos and Tsoubelis [25| we can gener- 
alize the energy conditions for a viscous anisotropic fluid. 

For the energy- momentum tensor Segre type [111, 1] and if Ao denotes the eigenvalue 
corresponding to the timelike eigenvector, the general energy conditions are equivalent to 
the following relations between the eigenvalues of the energy-momentum tensor: 



a) weak energy condition 



and 



b) dominant energy condition 



c) strong energy condition 



and 



A > 0, (68) 



- A + A, > 0, (69) 



Ao < Ai < -Ao, (70) 



A + E A *>°> ( 71 ) 



- Ao + A, > 0, (72) 

where the values % = 1,2,3 represent the eigenvalues corresponding to the spacelike eigen- 
vectors. 

The eigenvalues A of the energy-momentum tensor are the roots of the equation 

\T afj - \g a p\ = 0. (73) 
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Thus, we can rewrite equation ( 1731) as 



A 2 (fx + \) —ABq 
-ABq B 2 (p-\-(&)) 

o o c\ Vt - a - ce) 



o 
o 





C 2 (p t -\-(Q) 

where q = qB and the determinant of this equation is given by 

Xfi + \)(\-p + (Q) + q 2 ' x 
(\-p t + (Q) 2 A 2 B 2 C i = 0. 

Thus, one of the solutions of the equation (174"j) is 



(/i + A)(A-p + Ce) + g 



-2 



which can be rewritten as 



A 2 + (/i-P + ce)A + g 2 -Mp-ce) = o. 



The two roots of the equation (175]) are 



and 



A = --Gu-p + C© + A), 



A^-ifri-p + ce-A), 



(74) 



(75) 



(76) 



(77) 



(78) 



where 

A 2 = (^ + j9-Ce) 2 -4g 2 >0, (79) 

must be greater or equal to zero in order to have real solutions. This equation can be 
rewritten as 

\fM + p-CG\ -2\q\ >0. (80) 
The second solution of the equation (1741) is 



whose roots are given by 



(A- P , + ce) 2 = o, 



A 2 = A 3 = p t - C@- 



(81) 



(82) 
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A. Weak Energy Conditions 

From equations (168)) and (1771) we get the first weak energy condition written as 

/i-p + C© + A>0. (83) 

From equation (169|) . setting i = 1 and using equations ( 177)) and (ITS]) we get the second 
weak energy condition given by 

A > 0, (84) 



which is equal to the condition fl79|) . 

From equation (1691) . now setting * = 2, 3 (since A 2 = A 3 ) and using equations (!77|) and 
( I82p we get the third weak energy condition given by 

H - v + 2pt - C© + A > 0. (85) 
B. Dominant Energy Conditions 

From equation (1701) . setting i = 1 and using equations (177)) and ( 175]) we get the inequality 

-(^-p + C© + a) < -(^-p + ce-A) </i- P + ce + A, (86) 

which can be split into two inequalities, given by 

A > 0, (87) 

and 

fi-p + (Q>0. (88) 

From equation (1701) . setting i = 2, 3 and using equations (177)) and (1521) we get the in- 
equality 

- in - P + C© + A) < 2 (ft - C©) < /i - P + C© + A, (89) 
which again we can split it into two inequalities, given by 

li - p + 2p t - C© + A > 0, (90) 

and 

/i-p-2p t + 3C© + A>0. (91) 
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C. Strong Energy Conditions 



Substituting equations ( 1771) . (ITS]) and (182]) into equation (ITT]) we get the first strong energy 
condition given by 

2 Pt - 2(6 + A > 0. (92) 

Since one of the weak energy conditions, equation (1691) . is the same for the strong energy 
condition [equation (175]) ]. thus we have that the second and third strong energy conditions 
are equal to equations fl84l - (l85]) . given by 

A > 0, (93) 

and 

fj, - p + 2p t - (6 + A > 0. (94) 
D. Summary of the Energy Conditions 

Summarizing the results, we rewrite the energy conditions. The energy conditions for 
a spherically symmetric fluid whose energy-momentum tensor is given by equation ([3]) are 
fulfilled if the following inequalities are satisfied: 

(z) |/i + p-Ce|-2]g|>0, (95) 
(ii) [i - p + 2p t + A - (6 > 0, (96) 

and besides, 

a) for the weak energy conditions 

(Hi) /i-p + A + C6>0, (97) 

b) for the dominant energy conditions 

(iv) /U-p + C@>0, (98) 
(v) ii - p - 2p t + A + 3(6 > 0, (99) 
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c) for the strong energy conditions 



(vi) 2 Pt + A- 2(G > 0, (100) 

where A = JJjj, + p — C@) 2 — 4g 2 . 

In order to verify the energy conditions, we have plotted the time evolution of all the 
conditions, for several radii and for two values of ( (0 and 1000), as we can see in the figures 
IU [10] and [UJ For the sake of comparison with the model ( ^ 0, we have plotted all the 
conditions (l95|)- ffT00|) for C = 0. 

From the figures M^i) andfTOTi) we can conclude that only the inequality +p — (Q\ — 
2\q\ > 0] is not satisfied during all the collapse and for any radius. This inequality is not 
satisfied for the innermost radii (r < 0.2rs) and for the latest stages of the collapse. The 
condition (11001) is not satisfied for r < 0.2r s [figure [TTlfai)] because the inequality (|95|) 
[A > 0] is not satisfied for these radii and for the latest stages of the collapse. 



VII. PHYSICAL RESULTS 

As in previous papers lOl . [12|, we have calculated several physical quantities, 

as the total energy entrapped inside the S surface, the total luminosity perceived by an 
observer at rest at infinity and the effective adiabatic index, and we have compared them to 
the respective non-viscous ones. 

From equation fl35l) we can write using fj44|) - fr46|) and fl6T|) - fl66|) that 



m 
m 



f 



16r|(l -rl)h? 

A \2 



^(1 + ^)7^+ 



+44(1 - rlYih 2 ~ n + 16r|(l 



i) 2 / 2 



(101) 



where 



m 



r 2 B' + °- 

2B n 



(102) 



We can see from figure [12] that the total energy entrapped inside the hypersurface S 
vanishes at the time —2.0 x 10 -5 s, approximately. This means that the star radiates all its 
mass during the collapse and this explains why the apparent horizon never forms. We can 
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FIG. 9: The energy conditions (|95p - (|98p . for the model without bulk viscosity, where £ = 0. The 
time is in units of seconds and all the others quantities are in units of sec -2 . 
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FIG. 10: The energy conditions (|95p - (|98p . for the model without bulk viscosity, where £ = 1000. 
The time is in units of seconds and all the others quantities are in units of sec -2 . 
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FIG. 11: The energy conditions (|99p - (|100p . for the model with or without bulk viscosity, where 
C = and C = 1000. The time is in units of seconds and all the others quantities are in units of 
sec . 
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also observe from figure [121 that the mass inside £ is equal for both models, with or without 
bulk viscosity. This means that they radiate the same amount of mass during the evolution. 
In the figure [T51 we can see the evolution of the mass from a previous model [10|. In contrast 
of the result of this work, the former model radiates about 33% of the total mass of the star, 
before the formation of the black hole. 

Using the equations (|4T|) and (H4|) - fH6|) we can write the luminosity of the star as 



™§(i + 4) 4 / a „ 

K 3 X 



X 



l + rl 



A ' --riff 



2r|(l 



-C 



' 1 + 4 ' 
l-r|. 



' f h s 



(103) 



We can see from figure [TH that the luminosity perceived by an observer at rest at infinity 
increases exponentially until the time —2.0 x 10 -5 , when the total mass of the star vanishes. 
In the figure [TS] we can see the evolution of the luminosity from a previous model jlQ |. In 
contrast of the result of this work, the luminosity of the star also increases exponentially, 
reaching a maximum and after it decreases until the formation of the black hole. 

The effective adiabatic index can be calculated using the equations (l48]) - (pf9]) . (150"]) and 
fl6T]) - fl66|) . Thus, we can write that 



eff 



<9(lnp) 
<9(ln /i) 



- r=const 



P 



/i 



{{ [288r 2 e(r) + 12d{r)k{r)] f 2 + 6ac(r)k(r)hf f - 6K(A c(r)k(r)(h 3 + h)f 2 /a ] fh+ 



+ 



+c (r)k(r)j(r,t)f[3h 2 f + bh 2 (l-f) 
+k{r)hf \c{r){12bh 2 + aj{r,t)f 2 ) - 12d(r)h 2 ]} x 
X2" 1 {{24r 2 d{r)f 2 + k{r) [c{r){3h 2 f 2 + afhf + bh 2 {l - f 2 ))- 
- 2d{r)f 2 ] } fh + c{r)k{r){3h 3 f 3 + ah 2 f f 2 ) + 
+k{r)f [bc{r)h 3 {l - f 2 ) + 2d{r)h 3 ] - 2c{r)k{r)f 2 fhh{l + /i 2 )/a } _1 x 



x {l2r 2 d(r)f 2 + k(r) [c(r)h 2 f 2 + d(r)(h 2 - f 2 ) 
+ 2c(r)k(r)ffhh] x 



+ 
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-0.001 -0.0008 -0.0006 -0.0004 -0.0002 



z 

FIG. 12: Time behavior of the total energy entrapped inside the surface S for the models with or 
without bulk viscosity. The time, m and mo are in units of seconds. The symbols ( > mean that 
the plotted quantity is independent of (. 
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FIG. 13: Time behavior of the total energy entrapped inside the surface £ for the models with or 



without bulk viscosity, from a previous model 



10l | . The time, to and too are in units of seconds. 



The symbols ( > mean that the plotted quantity is independent of (. 
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FIG. 14: Time behavior of the luminosity perceived by an observer at rest at infinity for the models 
with or without bulk viscosity. The time is in units of second and the luminosity is dimensionless. 
The symbols ( > mean that the plotted quantity is independent of (. 
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FIG. 15: Time behavior of the luminosity perceived by an observer at rest at infinity for the models 
with or without bulk viscosity, from a previous model The time is in units of second and the 
luminosity is dimensionless. The symbols £ > mean that the plotted quantity is independent of 

C- 
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where 



x \l2r 2 e{r)f + jfe(r) {c(r)j(r, t)hff + 3[c(r)6 - d(r)} x 

x (h 2 -f 2 )}-c(r)k(r)l(r,t)f 2 h}'\ (104) 

c(r) =r 2 m 2 (l + r|) 8 , (105) 

d(r)=r|s 2 (r), (106) 

e(r) = r|(r|-r 2 )(7(r), (107) 

k(r) = (1 + r 2 ) 2 , (108) 

j(r,t) = 3a + 6/tCA /i. (109) 



and 



Z(r,t) = 3a o b-3n(A h. (110) 

Comparing the figures for r e g (£ = and ( ^ 0) we can see that the time evolution of 
the effective adiabatic indices are not very different graphically. This is the reason to plot 
the quantity 8T = T cS (( = 0) — T cS (( ^ 0) instead of T e g for the ( ^ models. We can 
note in figure [16] (( = 0) that shortly before the peak of luminosity (see figure fT4"j) there is 
a large discontinuity in r e fr due mainly to the behavior of the pressure. The effect of the 
viscosity is to increase much more these discontinuities. 

In the figure [T7| we can see the evolution of the effective adiabatic indices from a previous 
model [h3| . We can note comparing it with the figure [16] that the effective adiabatic index 
diminishes due to the bulk viscosity, thus increasing the instability of the system, in both 
models, in the former paper 10] and in this work. This characteristic might be model 
independent. 

Finally, models of radiating viscous spheres have been presented by Herrera, Jimenez 
and Barreto [26j . This work is particularly relevant for the proposed discussion because the 
conclusion concerning the effective adiabatic index is the same in both cases. Namely, an 
increasing of the critical adiabatic index required for stability [27], or equivalently, a decreas- 
ing of the effective adiabatic index, induced by viscosity. Since the models considered in each 
case are completely different, we suggest that this effect seems to be model independent. 



34 



100 



50 



i i i i | i i i r~ 

t = (This work) 

r = 10- 10 

i = 0.2 i \ 

r = 0.4 r 1 

r = 0.6 r l 

r = 0.8 r. 







-50 



-100 



(- = 500 (This work) 

r = 10- 10 r z 

r = 0.2 i y 

0.4 r v 



r 
r 
r 



0.6 r l 
0.8 r. 



I 1 



: 







-0.003 



i i i i | i i i i |~ 

t = 100 (This work) 
r = 10- 10 r v 



0.2 r z 

0.4 r Y 

0.6 r x 

r = 0.8 r v 



r 
r 
r 



t = 1000 (This work) 
r = 10- 10 r. 



r = 0.2 r E 

r = 0.4 r E I 

r = 0.6 r E / 

r = 0.8 r E / 







I I I I I I I I I I I I' I EE I I I I I I I I I I I I i QL . II I 







-0.002 -0.001 
t 







-0.002 -0.001 
t 







FIG. 16: Time behavior of the effective adiabatic index r e fj for four values of £. The quantity 5T is 
defined as r c fj(£ = 0) — T e g(( / 0). The time is in units of seconds, T e fj and 8F are dimensionless. 
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FIG. 17: Time behavior of the effective adiabatic index T e ff for four values of from a previous 
model [l(J. The quantity 5T is defined as T c g (( = 0) — r c fj (( ^ 0). The time is in units of seconds, 
r e fj and 5T are dimensionless. 
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VIII. CONCLUSIONS 



A new model is proposed to a collapsing star consisting of an anisotropic fluid with 



bulk viscosity, radial heat flow and outgoing radiation. In a previous paper [11] one of us 
has introduced a function time dependent into the g rr , besides the time dependent metric 
functions gee and g^. We have generalized this previous model by introducing bulk viscosity 
and we have compared it to the non- viscous collapse. 

The behavior of the density, pressure, mass, luminosity and the effective adiabatic index 
was analyzed. We have also compared to the case of a collapsing fluid with bulk viscosity 
of another previous model [10], for a star with 6 M . 

As we have shown the black hole is never formed because the apparent horizon formation 
condition is never satisfied. This could be interpreted as the formation of a naked singularity, 
as Joshi, Dadhich and Maartens jl4] have suggested. However this is not the case because 
the star radiates all its mass before it reaches the singularity at r = and t = 0. Not even 
a marginally naked singularity is formed by the same reason, since in this case the apparent 
horizon should coincide with the singularity at r = and t = 0. 

The density and pressure have negative values although physically this could be consid- 
ered unreasonable. However, due to the heat flow (the term A) the energy conditions are 
partially satisfied. 

The pressure of the star, at the beginning of the collapse, is isotropic but due to the 
presence of the bulk viscosity the pressure becomes more and more anisotropic. 

The star radiates all its mass during the collapse and this explains why the apparent 
horizon never forms. In contrast of the result of this work, the former model radiates about 
33% of the total mass of the star, before the formation of the black hole. 

An observer at infinity will see a radial point source radiating exponentially until reaches 
the time of maximum luminosity and suddenly the star turns off because there is no more 
mass in order to be radiated. In contrast of the former model 1(| where the luminosity also 
increases exponentially, reaching a maximum and after it decreases until the formation of 
the black hole. 

The effective adiabatic index has a very unusual behavior because we have a non-adiabatic 
regime in the fluid due to the heat flow. The index becomes negative since the hydrodynamic 
pressure and the density may become negative. Besides, in this case, neither the density is 
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the measure of the total energy densit y o f a given piece of matter nor the hydrodynamic 



pressure the only opposing contraction 
the bu 
paper 



j. The effective adiabatic index diminishes due to 
k viscosity, thus increasing the instability of the system, in both models, in the former 
lo| and in this work, showing that this characteristic might be model independent. 
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